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Abstract 

We  present  an  algorithm  for  locating  and  orienting  a  set  of  sensor  arrays  that  have  been  deployed 
in  a  scene  at  unknown  locations  and  orientation  angles.  This  self-calibration  problem  is  solved 
using  a  number  of  source  signals  also  deployed  in  the  scene.  We  assume  each  array  can  estimate 
the  time-of-arrival  and  direction-of-arrival  (with  respect  to  the  array’s  local  orientation  coordi¬ 
nates)  of  every  source.  From  this  information  we  compute  the  array  locations  and  orientations. 
We  consider  four  subproblems,  in  which  the  source  signals  or  emission  times  are  either  known  or 
unknown.  We  develop  necessary  conditions  for  solving  the  self-calibration  problem  and  provide 
a  maximum  likelihood  solution  and  corresponding  location  error  estimate. 


1.  Introduction 


Unattended  Ground  Sensors  (UGSs)  arc  becoming  increasingly  important  for  providing  situational  awareness 
in  battlefield  deployments  [  1  ] .  The  basic  concept  is  to  deploy  a  large  number  of  low-cost,  self-powered  sensors 
that  acquire  and  process  data.  The  sensors  typically  consist  of  an  array  of  microphones  to  detect,  track,  and 
classify  acoustic  signatures.  In  addition,  seismic  and  low-cost  imaging  sensors  may  also  be  present. 

We  consider  a  sensor  deployment  architecture  as  shown  in  Figure  1 .  A  number  of  low-cost  sensors,  each 
equipped  with  a  local  processor,  a  low-power  communication  transceiver,  and  one  or  more  sensing  capabil¬ 
ities,  is  set  out  in  a  region.  Sensor  elements  may  collect  acoustic,  seismic,  and/or  image  data.  Each  sensor 
monitors  its  environment  to  detect,  track,  and  characterize  signatures.  The  sensed  data  is  processed  locally, 
and  the  result  is  transmitted  to  a  local  Central  Information  Processor  (CIP)  through  a  low-power  communi¬ 
cation  network.  The  CIP  fuses  sensor  information  and  transmits  the  processed  information  to  a  more  distant 
command  center. 

In  order  to  fuse  sensor  information  at  the  CIP  or  command  center,  it  is  important  to  know  the  location  and 
orientation  of  each  sensor.  Ground  sensors  arc  placed  in  the  field  by  persons,  by  an  air  drop,  or  by  artillery 
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Figure  1 :  Sensor  network  architecture.  A  number  of  low-cost  sensors  arc  deployed  in  a  region.  Each  sensor 
communicates  to  a  local  CIP,  which  relays  information  to  a  more  distant  command  center. 


launch.  For  careful  hand  placement,  accurate  location  and  orientation  of  the  sensors  can  be  assumed; 
however,  for  many  sensor  deployment  situations,  it  is  difficult  or  impossible  to  know  accurately  the  location 
and  orientation  of  each  sensor.  One  could  equip  each  sensor  with  a  GPS  and  compass  to  obtain  location  and 
orientation  information,  but  this  adds  to  the  expense  and  power  requirements  of  the  sensor  and  may  increase 
susceptibility  to  jamming.  Thus,  there  is  interest  in  developing  methods  to  self-calibrate  the  sensor  array 
with  a  minimum  of  additional  hardware  or  processing. 

In  this  paper  we  consider  an  approach  to  array  self-calibration  using  sources  in  the  field.  A  number  of 
signal  sources  arc  deployed  in  the  same  region  as  the  sensors  (see  Figure  2).  Each  source  generates  a  unique 
signature  that  is  detected  by  the  sensors.  From  the  time-of-arrival  (TOA)  and  direction-of-arrival  (DOA)  of 
each  source  signal,  we  compute  the  unknown  locations  and  orientations  of  the  sensors.  We  consider  four 
related  subproblems,  in  which: 

1.  the  source  locations  and  emission  times  are  known, 

2.  the  source  locations  arc  known  and  emission  times  arc  unknown, 

3.  the  source  locations  arc  unknown  and  emission  times  arc  known, 

4.  the  source  locations  and  emission  times  are  unknown. 

Several  researchers  have  considered  the  problem  of  array  calibration,  but  less  work  is  devoted  to  calibrating 
networks  of  sensors.  A  number  of  papers  have  considered  calibration  of  both  narrowband  and  broadband 
arrays  of  sensors  to  improve  direction-of-arrival  estimation  accuracy  [2,  3,  4,  5,  6,  7],  These  papers  assume 
knowledge  of  the  nominal  sensor  locations  high  (or  perfect)  signal  coherence  between  the  sensors.  Calibration 
requirements  for  acoustic  ground  sensors  arc  discussed  in  [8] .  Research  on  blind  beamforming  considers  a 
related  problem  of  forming  a  maximum  power  beam  to  a  source  without  computing  the  source  locations  [9], 
A  recent  paper  considers  sensor  self-calibration  using  a  single  acoustic  source  that  travels  in  a  straight  line 
[10]. 

An  outline  of  the  paper  is  as  follows.  In  Section  2  we  present  a  statement  of  the  problem  and  justify  our 
assumptions.  In  Section  3  we  first  consider  necessary  conditions  for  a  self-calibration  solution  and  present 
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Figure  2:  Array  self-calibration  scenario. 


methods  for  solving  the  self-calibration  problem  with  a  minimum  number  of  sensors  and  sources.  These 
methods  provide  initial  estimates  for  an  iterative  descent  computation  needed  to  obtain  maximum  likelihood 
calibration  parameter  estimates  derived  in  Section  4.  Bounds  on  the  calibration  uncertainty  arc  also  derived. 
Section  5  presents  numerical  examples  to  illustrate  the  approach,  and  Section  6  presents  conclusions  and 
directions  for  future  work. 


2.  The  Self-Calibration  Problem 


Assume  we  have  a  set  of  A  sensors,  each  with  unknown  location  {a*  =  ( xr ,  and  unknown  orientation 

angle  0;  with  respect  to  a  reference  direction  ( e.g .,  North).  We  consider  the  two-dimensional  problem  in 
which  the  sensors  he  in  a  plane  and  the  unknown  reference  direction  is  azimuth;  an  extension  to  the  three- 
dimensional  case  is  possible  using  similar  techniques. 

In  the  array  held  arc  also  placed  S  point  source  signals  at  locations  {sj  =  ( Xj,yj)}j=1 .  The  source  locations 
may  be  known  or  unknown.  Each  source  emits  a  finite-length  signal  that  begins  at  time  tj ;  the  emission  times 
may  be  known  or  unknown.  We  thus  consider  four  related  subproblems,  depending  on  the  prior  knowledge 
of  the  source  locations  and  their  emission  times. 

We  initially  assume  each  emitted  source  signal  is  detected  by  all  of  the  sensors  in  the  held  (the  extension  to 
partial  measurements  is  considered  in  Section  4.4)  and  that  each  sensor  measures  the  time-of-arrival  (TOA) 
and  direction-of-arrival  (DOA)  for  that  source.  We  denote  the  measured  TOA  of  source  j  to  sensor  i  as  tr] 
and  the  measured  DOA  as  Oij. 

The  DOA  measurements  are  made  with  respect  to  a  local  (to  the  sensor)  frame  of  reference.  The  absolute 
directions  of  arrival  are  not  available  because  the  orientation  angle  of  each  array  is  unknown  (and  must  be 
estimated  in  the  calibration  procedure).  The  times  of  arrival  arc  measured  with  respect  to  a  known,  common 
time  base.  The  time  base  can  be  established  either  by  using  the  electronic  communication  network  linking 
the  sensors  or  by  synchronizing  the  sensor  processor  clocks  before  deployment.  The  time  base  needs  to  be 
accurate  to  a  number  on  the  order  of  the  sensor’s  time  of  arrival  measurement  uncertainty  (1  msec  in  the 


examples  in  Section  5). 

The  set  of  2 ,4 ,5’  measurements  arc  gathered  in  a  vector 


A'  = 


vec(T) 

vec(0) 


(2AS  x  1) 


where  vec(M)  stacks  the  elements  of  a  matrix  M  columnwise  and  where 
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Each  array  transmits  its  25  TOA  and  DOA  measurements  to  a  central  information  processor,  and  these  2  AS 
measurements  form  the  data  with  which  the  CIP  computes  the  array  calibration.  Note  that  the  communication 
cost  to  the  CIP  is  low,  and  the  calibration  processing  is  entirely  performed  by  the  CIP. 


Define  the  parameter  vectors 


[/?r,7T 
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We  denote  the  actual  TOA  and  DOA  of  source  signal  j  at  sensor  i  as  nj(a)  and  4>ij(a),  respectively,  and 
include  their  dependence  on  the  parameter  vector  a;  they  arc  computed  as: 


Tij(a)  =  tj  +  \\cii  -  Sj\\/c  (6) 

=  6i  +  Z(aj,  Sj)  (7) 

where  ||  •  ||  is  the  Euclidean  norm,  Z(£,  //)  is  the  angle  between  the  points  £,  ij  €  'R1.  and  c  is  the  signal 
propagation  velocity. 

Each  element  of  X  has  measurement  uncertainty;  we  model  the  uncertainty  as 

X  =  n(a)  +  E  (8) 

where  n{a)  is  the  noiseless  measurement  vector  whose  elements  arc  given  by  equations  (6)  and  (7)  for  values 
of  i.  j  that  correspond  to  the  vector  stacking  operation  in  (1),  and  where  E  is  a  random  vector  with  known 
probability  density  function. 

The  self-calibration  problem,  then,  is  given  the  measurement  X,  estimate  ft.  Note  that  none,  some,  or  all  of 
the  parameters  in  7  may  be  known  depending  on  the  particular'  subproblem  of  interest. 


3.  Existence  and  Uniqueness  of  Solutions 


In  this  section  we  address  the  existence  and  uniqueness  of  solutions  to  the  self-calibration  problem  and 
establish  the  minimum  number  of  sensors  and  sources  needed  to  obtain  a  solution.  We  assume  that  every 
sensor  detects  every  source  and  measures  both  TOA  and  DOA.  In  addition,  we  assume  the  TOA  and  DOA 
measurements  arc  noiseless.  We  establish  the  minimum  number  of  sources  and  sensors  needed  to  compute 
a  unique  calibration  solution  and  give  algorithms  for  finding  the  self-calibration  solution  in  the  minimal 
cases.  These  algorithms  provide  initial  estimates  to  an  iterative  descent  algorithm  for  the  practical  case  of 
non-minimal,  noisy  measurements  presented  in  Section  4. 

The  four  cases  below  make  different  assumptions  about  the  source  signal  parameters.  In  all  four  cases  the 
number  of  measurements  is  2AS ,  and  determination  of  ft  involves  solving  a  nonlinear  set  of  equations  for 
its  3 A  unknowns.  Depending  on  the  case  considered,  we  may  also  need  to  estimate  the  unknown  nuisance 
parameters  in  7.  The  result  in  each  case  is  summarized  in  Table  1 . 


Case  1:  Known  source  locations  and  emission  times. 

A  unique  solution  for  ft  can  be  found  for  any  number  of  sensors  as  long  as  there  arc  S  >  2  sources.  In  fact,  the 
location  and  orientation  of  each  sensor  can  be  computed  independently  of  other  sensor  measurements.  The 
location  of  the  ith  sensor,  cq,  is  found  from  the  intersection  of  two  circles  with  centers  at  the  source  locations 
and  with  radii  ( tn  —  tft)/c  and  ( ta  —  t\)/c.  The  intersection  is  in  general  two  points;  the  correct  location  can 
be  found  using  the  sign  of  622  —  On.  We  note  that  the  two  circle  intersections  can  be  computed  in  closed-form. 
From  the  known  source  and  sensor  locations  and  the  DOA  measurements,  the  sensor  orientations  can  also 
be  uniquely  found. 


Case  2:  Known  source  locations,  unknown  emission  times. 

For  5  >  3  sources  the  location  and  orientation  of  each  sensor  can  be  computed  in  closed  form  independently 
of  other  sensors.  A  solution  procedure  is  as  follows.  Consider  the  pair  of  sources  (si,  S2).  Sensor  i  knows 
the  angle  Oa  —  On  between  these  two  sources.  The  set  of  all  possible  locations  for  sensor  i  is  an  arc  of  a 
circle  whose  center  and  radius  can  be  computed  from  the  source  locations  (see  Figure  3).  Similarly,  a  second 
circular  arc  is  obtained  from  the  source  pair  (si,  S3).  The  intersection  of  these  two  arcs  is  a  unique  point  and 
can  be  computed  in  closed  form.  Once  the  sensor  location  is  known,  its  orientation  is  readily  computed  from 
any  of  the  three  DOA  measurements. 

A  solution  for  Case  2  can  also  be  found  using  5  =  2  sources  and  A  =  2  sensors.  The  solution  requires 
a  one-dimensional  search  of  a  parameter  over  an  finite  interval.  The  known  location  of  si  and  S2  and  the 
known  angle  On  —  On  means  that  array  01  must  lie  on  a  known  circular  arc  as  in  Figure  3.  Each  location 
along  the  arc  determines  both  source  emission  times.  At  exactly  one  position  along  the  arc,  the  emission 
times  are  consistent  with  the  measurements  from  the  second  sensor. 


Table  1:  Minimal  Solutions  for  Array  Self-Calibration 


Case 

#  Unknowns 

Minimum  A ,  5 

Comments 

Known  Locations 

Known  Times 

3A 

A  =  l,  S  =  2 

closed-form  solution 

Known  Locations 

Unknown  Times 

3  A  +  S 

A  =  l,  S  =  3 

closed-form  solution 

3  A  +  S 

A  =  2,  5  =  2 

1-D  iterative  solution 

Unknown  Locations 

Known  Times 

3{A-  1)  +2S 

A  =  2,  5  =  2 

closed-form  solution 

Unknown  Locations 
Unknown  Times 

3{A+S-  1) 

A  =  2,  5  =  3  or 

A  =  3,  5  =  2 

2-D  iterative  solution 

Case  3:  Unknown  source  locations,  known  emission  times. 

In  this  case  and  in  Case  4  below,  the  calibration  problem  can  only  be  solved  to  within  an  unknown  translation 
and  rotation  of  the  entire  sensor-source  scene  because  any  translation  or  rotation  of  the  entire  scene  does  not 
change  the  tij  and  9ij  measurements.  To  eliminate  this  ambiguity,  we  assume  the  location  and  orientation  of 
the  first  sensor  arc  known;  without  loss  of  generality  we  set  x\  =  y\  =8 1  =  0.  We  solve  for  the  remaining 
3(A  —  1)  parameters  in  (3. 

For  the  case  of  unknown  source  locations,  a  unique  solution  for  (3  is  computable  in  closed  form  for  5  =  2 
and  any  A  >  2  (the  case  A  =  1  is  trivial).  From  sensor  ai  the  range  to  each  source  can  be  computed  from 
Tj  =  (tij  —  tj)/c,  and  its  bearing  is  known,  so  the  locations  of  the  two  sources  can  be  found.  The  locations 
and  orientations  of  the  remaining  sensors  is  then  computed  using  the  method  of  Case  1. 


Case  4:  Unknown  source  locations  and  emission  times. 

For  this  case  it  can  be  shown  that  an  infinite  number  of  calibration  solutions  exists  for  A  =  5  =  2,1  but  that 
a  unique  solution  exists  in  almost  all  cases  for  either  A  =  2,  S  =  3  or  ,4  =  3,  S  =  2.  In  some  degenerate 
cases,  not  all  of  the  7  parameters  can  be  uniquely  determined,  although  we  do  not  know  of  a  case  for  which 
the  f3  parameters  cannot  be  uniquely  found. 

Closed-form  calibration  solutions  arc  not  known  for  this  case,  but  solutions  that  require  a  two-dimensional 
search  can  be  found.  We  outline  one  such  solution  that  works  for  either  A  =  2  and  5  >  3  or  5  =  2  and 
A  >  3.  Assume  as  before  that  sensor  ai  is  at  location  (x\ ,  y\ )  =  (0,0)  with  orientation  9\  =  0.  If  we 
knew  the  two  source  emission  times  t\  and  t,2,  we  can  find  the  locations  of  sources  .S|  and  S2  as  in  Case  3. 
From  the  two  known  source  locations,  all  remaining  sensor  locations  and  orientations  can  be  found  using 
the  procedure  in  Case  1,  and  then  all  remaining  source  locations  can  be  found  using  triangulation  from  the 
known  arrival  angles  and  known  sensor  locations.  The  solutions  will  be  inconsistent  except  for  the  correct 
values  of  t\  and  to-  The  calibration  procedure,  then,  is  to  iteratively  adjust  t\  and  to  to  minimize  the  error 

1Note  that  for  A  =  S  —  2  there  are  8  measurements  and  9  unknown  parameters.  The  set  of  possible  solutions  in  general  lies  on  a 
one-dimensional  manifold  in  the  9-dimensional  parameter  space. 


S1  s2 


Figure  3:  A  circular  arc  is  the  locus  of  possible  sensor  locations  whose  angle  between  two  known  points  is 
constant. 

between  computed  and  measured  time  delays  and  arrival  angles. 


4.  Maximum  Likelihood  Self-Calibration 


In  this  section  we  derive  a  maximum  likelihood  (ML)  estimator  for  the  unknown  array  location  and  orientation 
parameters.  The  algorithm  involves  the  solution  of  a  set  of  nonlinear  equations  for  the  unknown  parameters 
(and  the  unknown  nuisance  parameters  in  7).  The  solution  is  found  by  iterative  minimization  of  a  cost 
function;  we  use  the  methods  in  Section  3  to  initialize  the  iterative  descent.  In  addition,  we  derive  the 
Cramer-Rao  Bound  (CRB)  for  the  variance  of  the  unknown  parameters  in  a;  the  CRB  also  gives  high-SNR 
parameter  variance  of  the  ML  parameter  estimates. 


4. 1  The  Maximum  Likelihood  Estimate 


We  assume  the  measurement  uncertainty  E  in  equation  (8)  is  Gaussian  with  zero  mean  and  known  covariance 
E.  In  this  case  the  likelihood  function  is 


/(X;a)  =  (2^7^expH0<Xia)}  (9) 

Q(X;a )  =  [X  -  ^(a)]rE"1[X  - //(a)]  (10) 


A  special  case  is  when  the  measurement  errors  are  uncorrelated  and  the  TOA  and  DOA  measurement  errors 
are  Gaussian  with  zero  mean  and  variances  of  and  of,  respectively;  then  equation  (10)  becomes 
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In  the  four  cases  considered  in  Section  3,  some  of  the  parameters  in  a  arc  known.  We  denote  07  to  be  the 
unknown  parameters  in  a  and  cm  to  be  the  known  parameters  for  the  particular  case  of  interest.  Using  this 


notation  along  with  equation  (9),  the  maximum  likelihood  estimate  of  a\  is 


&i ml  =  argmax/(X,  «)  =  argmin  Q(X]  a) 

a.  i  a  i 


(12) 


4.2  Nonlinear  Least  Squares  Solution 

The  solution  of  (12)  involves  solving  a  nonlinear  least  squares  problem.  A  standard  iterative  descent  procedure 
can  be  used,  initialized  using  one  of  the  solutions  in  Section  3.  In  our  implementation  we  used  the  Matlab 
function  lsqnonlin. 

The  straightforward  nonlinear  least  squares  solution  we  adopted  converged  quickly  (in  several  seconds  for 
all  examples  tested)  and  displayed  no  symptoms  of  numerical  instability;  however,  alternative  methods  for 
solving  equation  (12)  may  reduce  computation.  One  can  divide  the  parameter  set  and  iterate  first  on  the  sensor 
location  parameters  and  second  on  the  remaining  parameters.  Although  the  sensor  orientations  and  source 
parameters  depend  nonlinearly  on  the  sensor  locations,  computationally  efficient  approximations  exist  (see, 
e.g. ,  [1 1]),  so  the  computational  savings  of  lower  dimensional  searches  may  exceed  the  added  computational 
cost  of  iterations  nested  in  iterations  if  the  methods  are  tuned  appropriately.  Similarly,  one  can  view  the 
source  parameters  as  nuisance  parameters  and  employ  estimate-maximize  (EM)  algorithms  to  obtain  the  ML 
solution  [12], 


4.3  Estimation  Accuracy 

The  Cramer-Rao  Bound  (CRB)  gives  a  lower  bound  on  the  covariance  of  any  unbiased  estimate  of  a \ .  It  is 
a  tight  bound  in  the  sense  that  a\ tML  has  parameter  uncertainty  given  by  the  CRB  for  high  signal-to-noise 
ratio;  that  is,  as  max,;  —y  0.  Thus,  the  CRB  is  a  useful  tool  for  analyzing  calibration  uncertainty. 

The  CRB  can  be  computed  from  the  Fisher  Information  Matrix  of  a\.  The  Fisher  Information  Matrix  is 
given  by  [13] 

Iai  =  E  { [Vai  In  f(T,  0;  a)}  [Vai  In  f(T,  0;  a)]T} 

The  partial  derivatives  arc  readily  computed  from  equations  (9),  (6),  and  (7);  we  find  that 

Iai  =  [G"(ai)]TS-1[G'(ai)]  (13) 

where  G'(a\)  is  the  2 AS  x  dim(ai)  matrix  whose  ijth  element  is  d Hi{ai) / d(ot\) j . 

For  Cases  3  and  4,  the  Fisher  information  matrix  is  rank  deficient  due  to  the  translational  and  rotational 
ambiguity  in  the  self-calibration  solution  in  those  two  cases.  In  this  situation,  two  approaches  can  be  taken. 

First,  one  can  assume  some  of  the  sensor  parameters  arc  known.  Examples  arc  to  assume  known  location 
and  orientation  of  a  single  sensor,  or  to  assume  known  location  of  two  sensors.  These  assumptions  might  be 
realized  by  equipping  one  sensor  with  a  GPS  and  a  compass,  or  by  equipping  two  sensors  with  GPSs.  Let 
d\  denote  the  vector  obtained  by  removing  these  assumed  known  parameters  from  a  i .  To  compute  the  CRB 
matrix  for  di  in  this  case,  we  first  remove  all  rows  and  columns  in  Iai  that  correspond  to  the  assumed  known 
parameters  then  invert  the  remaining  matrix  [13]: 


(14) 


The  second  approach  is  to  compute  the  CRB  of  the  parameter  vector  a\  subject  to  knowledge  of  the  translation 
and  rotation.  To  do  so  we  compute  an  eigenvalue  decomposition  of  Iai : 
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Except  in  degenerate  cases  (in  which,  for  example,  all  sensors  and  at  least  one  source  signal  are  col  linear),  it 
can  be  shown  that  U2  has  3  columns  and  that  its  columns  span  the  subspace  corresponding  to  overall  scene 
translation  and  rotation.  Then  the  constrained  CRB  of  the  parameter  vector  a\  subject  to  knowledge  of  the 
translation  and  rotation  is  given  by  the  pseudoinverse  of  Iai  [14]: 

Ca  =  (/a)f  =  U^Uf  (16) 


4.4  Partial  Measurements 


So  far  we  have  assumed  that  every  sensor  detects  and  measures  both  the  TOA  and  DOA  from  every  source 
signal.  In  this  section  we  relax  that  assumption.  We  assume  each  emitted  source  signal  is  detected  by  only 
a  subset  of  the  sensors  in  the  field  and  that  a  sensor  that  detects  a  source  may  measure  the  time-of-arrival 
(TOA)  and/or  the  Direction-of-Arrival  (DOA)  for  that  source,  depending  on  its  capabilities.  We  denote  the 
availability  of  a  measurement  using  two  indicator  functions  /A  and  ifj,  where 

44  €{0,1}  (17) 

If  sensor  i  measures  the  TOA  (DOA)  for  source  j,  then  /A  =  1  (ifj  =  1);  otherwise,  the  indicator  function 
is  set  to  zero.  Furthermore,  let  L  denote  the  2 ,4, S'  x  1  vector  whose  fcth  element  is  1  if  A/,  is  measured  and  is 
0  if  Xj.  is  not  measured;  L  is  thus  obtained  by  forming  A  x  S  matrices  I1  and  Ie  and  stacking  their  column 
into  a  vector  as  in  equation  (1).  Finally,  define  X  to  be  the  vector  formed  from  elements  of  X  for  which 
measurements  arc  available,  so  Xj.  is  in  A  if  L\.  =  1. 


The  maximum  likelihood  estimator  for  the  partial  measurement  case  is  similar  to  equation  (12)  but  uses  only 
those  elements  of  A  for  which  the  corresponding  element  of  L  is  one.  Thus, 

cti, ml  =  arg  min  Q( A;  a)  (18) 


where  (assuming  uncorrelated  measurement  errors  as  in  equation  (11)), 

A  S 
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The  Fisher  Information  Matrix 
available  measurements;  thus 


for  this  case  is  similar  to  equation  (13),  but  includes  only  information  from 


Iai  =  [G"(ai)]TS-1[G"(ai)] 


(20) 
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where 


(21) 


We  note  that  when  partial  measurements  arc  available,  the  ML  calibration  may  not  be  unique.  For  example, 
if  only  TOA  measurements  arc  available,  a  scene  calibration  solution  and  its  mirror  image  have  the  same 
likelihoods.  A  complete  understanding  of  the  uniqueness  properties  of  solutions  in  the  partial  measurement 
case  is  a  topic  of  current  research. 


5.  Numerical  Results 


We  present  some  numerical  examples  of  the  self-calibration  procedure.  Ten  sensor  arrays  and  eleven  sources 
arc  randomly  placed  in  a  2  km x 2  km  region;  the  sensor  orientations  and  source  emission  times  arc  also 
randomly  chosen.  Figure  4  shows  the  locations  of  the  sensors  and  sources.  We  assume  every  sensor  detects 
each  source  emission  and  measures  the  TOA  and  DOA  of  the  source.  The  measurement  uncertainties  are 
Gaussian  with  standard  deviations  of  at  =  1  msec  for  the  TOAs  and  a$  =  3°  for  the  DOAs.  Neither  the 
locations  nor  emission  times  of  the  sources  arc  assumed  to  be  known,  and  to  eliminate  the  translation  and 
rotation  uncertainty  in  the  scene,  we  assume  either  two  sensors  have  known  locations  or  one  sensor  has 
known  location  and  orientation. 

Figure  4  also  shows  the  two  standard  deviation  (2a)  location  uncertainty  ellipses  for  both  the  sources  and 
sensors  assuming  the  locations  of  sensor  arrays  .41  and  ,42  arc  known.  These  ellipses  appeal-  as  small  dots  in 
the  figure;  an  enlarged  view  for  two  sensors  are  shown  in  Figure  5.  As  a  quantitative  measure  of  performance, 
we  compute  for  each  sensor  the  equivalent  2a  uncertainty  radius,  defined  as  the  geometric  mean  of  the  major 
and  minor  axis  lengths  of  the  2 a  uncertainty  ellipse.  We  find  that  the  average  of  the  2 a  uncertainty  radii  for 
all  ten  sensors  is  0.80  m  for  this  example. 

Figure  6  shows  a  similar  uncertainty  plot,  but  in  this  case  we  assume  that  both  the  location  and  orientation  of 
sensor  ,41  is  known.  In  comparison  with  Figure  4,  we  see  much  larger  uncertainty  ellipses  for  the  sensors, 
especially  in  the  direction  tangent  to  circles  with  center  at  sensor  Al.  The  average  2a  uncertainty  radius  is 
3.28  m  in  this  case.  The  high  tangential  uncertainty  is  primarily  due  to  the  DOA  measurement  uncertainty 
with  respect  to  a  known  orientation  of  sensor  Al.  By  comparing  Figures  4  and  6,  we  see  that  it  is  more 
desirable  to  know  the  locations  of  two  sensors  than  to  know  the  location  and  orientation  of  a  single  sensor; 
thus,  equipping  two  sensors  with  GPS  systems  results  in  lower  uncertainty  than  equipping  one  sensor  with 
a  GPS  and  a  compass. 

Figure  7  shows  the  effect  of  increasing  the  number  of  sources  on  the  average  2er  uncertainty  radius.  We  plot 
the  uncertainties  found  from  using  from  2  through  11  sources,  starting  initially  with  sources  51  and  52  and 
adding  sources  53,  54,  ...  511  at  each  step.  We  see  the  uncertainty  reduces  dramatically  when  the  number 
of  sources  increases  from  2  to  3  and  then  decreases  more  gradually  as  more  sources  are  added.  For  this 
example,  an  average  2cr  uncertainty  radius  of  1-3  meters  is  obtained  when  more  than  five  signal  sources  are 
used  for  calibration. 

The  algorithms  require  low  communication  overhead  as  each  sensor  needs  to  communicate  only  25  scalar 
values  to  the  CIP.  Computation  of  the  calibration  solution  takes  place  at  the  CIP;  for  the  examples  presented 
the  calibration  computation  takes  on  the  order  of  10  seconds  using  Matlab  on  a  standard  personal  computer. 


6.  Conclusions 


We  have  presented  a  procedure  for  calibrating  the  locations  and  orientations  of  a  network  of  sensors.  The 
calibration  procedure  uses  source  signals  that  are  placed  in  the  scene  and  computes  array  and  source  unknowns 
from  estimated  time-of-arrival  and/or  direction-of-arrival  estimates  obtained  for  each  source-sensor  pair.  We 
present  maximum  likelihood  solutions  to  four  variations  on  this  problem,  depending  on  whether  the  source 
locations  and  signal  emission  times  arc  known  or  unknown.  We  also  discuss  existence  and  uniqueness 
of  solutions  and  algorithms  for  initializing  the  nonlinear  minimization  step  in  the  maximum  likelihood 
estimation.  An  analytical  expression  for  the  sensor  location  and  orientation  error  covariance  matrix  is  also 
presented  for  each  of  the  four  problem  variations.  A  maximum  likelihood  calibration  algorithm  for  the  case 
of  partial  calibration  measurements  was  also  developed. 

The  algorithms  require  minimal  communications  from  the  sensors  to  a  Central  Information  Processor,  and 
computation  of  the  calibration  solution  takes  on  the  order  of  ten  seconds  using  Matlab  on  a  standard  personal 
computer  for  examples  considered. 
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Figure  4:  Example  scene  showing  ten  sensor  arrays  (stars)  and  eleven  sources  (squares).  Also  shown  are  the 
2ct  location  uncertainty  ellipses  of  the  sensors  and  sources;  these  are  on  average  less  than  1  m  in  radius  and 
show  as  small  dots.  The  locations  of  sensors  Al  and  .42  arc  assumed  to  be  known. 


Figure  5:  Two  standard  deviation  location  uncertainty  ellipses  for  sensors  A3  (left)  and  .49  (right)  from 
Figure  4. 
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Figure  6:  The  2cr  location  uncertainty  ellipses  for  the  scene  in  Figure  4  when  the  location  and  orientation  of 
sensor  Al  is  assumed  to  be  known. 


Figure  7:  Average  2a  location  uncertainty  radius  for  the  scenes  in  Figures  4  and  6  as  a  function  of  the 
number  of  source  signals  used. 


